Quantum criticality in disordered bosonic optical lattices 
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Using the exact Bose-Fermi mapping, we study universal properties of ground-state density dis- 
tributions and finite-temperature quantum critical behavior of one-dimensional hard-core bosons in 
trapped incommensurate optical lattices. Through the analysis of universal scaling relations in the 
quantum critical regime, we demonstrate that the superfluid to Bose glass transition and the general 
phase diagram of disordered hard-core bosons can be uniquely determined from finite-temperature 
density distributions of the trapped disordered system. 
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Introduction. - In recent years ultracold atomic systems 
have proven to be powerful quantum simulators for inves- 
tigating various challenging many-body physics. In com- 
parison with traditional condensed matter systems, both 
the potential and interaction of trapped atomic systems 
can be experimentally implemented in a controlled way, 
which leads to experimental breakthroughs in the study 
of fundamental model systems, including the realization 
of the Tonks-Girardeau (TG) gas and the observa- 
tion of Anderson localization in disordered bosonic op- 
tical lattices 0-0 • The interplay of disorder and inter- 
actions for bosons has often been studied in the context 
of Bose-Hubbard model with disordered or incommensu- 
rate potentials. Theoretical studies have indicated that 
a superfluid to Bose-glass (BG) transition occurs for in- 
creasing disorder [1, Despite intensive studies [8-ll|, 
no exact solutions for disordered interacting bosons are 
known except that in the TG limit the one-dimensional 
(ID) disordered Bose-Hubbard models are exactly solv- 
able via a Bose-Fermi mapping 12h14|. 

So far most of theoretical studies on ID disordered 
bosonic gases focused on zero temperature cases, as 
finite-temperature phase transitions are generally absent 
in ID systems. However, considering realistic experi- 
ments are always carried on at finite temperature for 
trapped systems, it is important to understand how to 
unambiguously determine the zero temperature phase di- 
agram from the knowledge of finite temperature density 
profiles of trapped gases [lUni. In the presence of a har- 
monic trap, Anderson plateaus are found in the average 
density profile of the hard-core bosons in the incommen- 
surate optical lattice [3]. The Anderson plateau can be 
viewed as the signature of Anderson localized states with 
superfluid to BG transitions taking place at edges of the 
plateau. At finite temperatures, Anderson plateaus are 
smeared out by thermal fluctuations and thus no sharp 
boundaries for superfluid and BG phases could be de- 
tected. In order to understand zero temperature phase 
transitions, we must consider the so-called quantum crit- 
ical regime where some universal scaling relations gov- 
ern the finite-temperature physics. In principle, one can 



extract the zero-temperature phase diagram of the dis- 
ordered Bose system from these universal relations. Un- 
fortunately, it is very hard to calculate these universal 
functions due to lack of reliable general techniques for dis- 
ordered many-body problems. Since quantum and ther- 
mal fluctuations strongly couple together in the quantum 
critical regime, highly precise calculation of the finite- 
temperature data is needed in order to extract correct 
zero-temperature phase diagram. In this work, by us- 
ing the exact numerical method based on the Bose-Fermi 
mapping, we can calculate finite-temperature properties 
of disordered hard-core bosons exactly, and thus demon- 
strate that the zero-temperature phase diagram can be 
uniquely determined from finite-temperature density dis- 
tributions of the trapped gas. 

To make progress in unveiling the quantum criticality 
in disordered Bose systems quantitatively, we study a 
disordered ID TG gas with the effect of disorder being 
mimicked by an incommensurate potential [H, 0] , which 
is described by 
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where b\ (bi) is the creation (annihilation) operator of the 
boson fulfilling the hard-core constraints, i.e., the on-site 
anticommutation ({hi, b\} — 1) and bj] = for i ^ j, 

and hi = SjSj. The hopping amplitude t is set to be the 
unit of the energy (t = 1), and V\ is given by 

Vi = Vicos{a2m + 8) + V H (i~i ) 2 - (2) 

Here Vi is the strength of incommensurate potentia with 
a being an irrational number and 5 an arbitrary phase 
which is chosen to be zero for convenience, and Vh is 
the strength of the additional harmonic trap with iq be- 
ing the position of trap center. In the TG limit, the 
Hamiltonian {T]) can be mapped to a noninteracting spin- 
less Fermi model, which allows us to calculate even the 
finite-temperature properties by using the exact numeri- 
cal method [3, EH. 
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FIG. 1: (Color online) Local average density profiles rela- 
tive to the scaled position for different systems with same 
characteristic density(p = 1). The three systems are: 200 
bosons, V H = 2.5 x 10" 5 ; 125 bosons, V H = 6.4 x 10" 5 ; and 
250 bosons, V H = 1.6 x 10~ 5 with 1000 sites, Vi = 1 and 
a = (\/E — l)/2. Insert: corresponding density profile for the 
system with 200 bosons. 
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FIG. 2: (Color online) Local average density profiles for sys- 
tems with different particle number (characteristic density). 
The other parameters of system are 1000 sites, Vi = 2.4, 
a = (-v/5- l)/2 and V H = 2.5 x 10" 5 . 



Universal properties of zero-temperature density 
distribution.- In order to characterize universal features 
of lattice systems in harmonic traps, we make use of 
the length scale £ = (Vh /t) -1 / 2 and the characteristic 
density p = N/( introduced in Ref.[20|. Without the 
incommensurate potential the system has been studied 
by Rigol and Muramatsu in detail [3], and they found 
that there is a critical characteristic density in the 
system (p c ~ 2.6 — 2.7). For low characteristic density( 
p < p c ), the whole system is in a uniform phase (super- 
fluid phase) at zero temperature, and for p > p c two 
different phases coexist with a Mott insulating plateau 
in the middle of the trap surrounded by superfluid 
phases on two sides. In the presence of incommensurate 
potential 14|, the density profiles basically still have 
the arc shape for weak Vr> but there are a lot of drastic 
oscillations induced by the incommensurate potential 
(see the insert picture in FigfT]). In order to reduce 
the drastic oscillations in density profiles, we define 
the local average density Hi = J2j=-M ni +j/(2M + 1), 
where 2M + 1 is the length to count the local average 
density with M < L and m = (*hcbI^I*hcb) witn 
^hcb * ne ground state of the Hamiltonian (|T|). When 
Vj —> 0, the local average density profile is almost the 
same as the density profile. As Vj increases, Anderson 
plateaus appear at shoulders of the arc, then become 
wider and wider but the height of the plateau does not 
change with V/. For sites in the Anderson plateau the 
one particle density matrices pij — (^hcbI^I^'I^hcb) 
have an exponential-law decay which is the character of 
the particle in the BG phase. While for sites outside 
the plateau, pij exhibit a power-law decay which is the 
character of the particle in the superfluid phase. 

Now we discuss the universal properties of Anderson 
plateaus. In FigfT] we show local average density pro- 



files for three systems with the same characteristic den- 
sity relative to the scaled position x = (i — io)/C- They 
are almost identical except of tiny differences caused by 
the setting of M = 10 for all three systems with differ- 
ent length scale £. So for systems with different par- 
ticle number and strength of the harmonic trap, An- 
derson plateaus are the same so long as these systems 
have the same characteristic density. From now on, we 
fix the strength of harmonic trap and change particle 
numbers to achieve different characteristic densities. By 
doing this we don't need to scale the position for com- 
parison among different systems. As shown in Fig [2] 
we plot local average density profiles for systems with 
fixed Vh and different particle numbers (different p). 
Here we choose Vj > 2 to generate as many Ander- 
son plateaus as possible in local average density pro- 
files. Due to the incommensurate potential taking the 
form of cos(27rai), for any a £ (—00,00) one can al- 
ways choose a corresponding a' £ [0, 0.5) which generates 
the same incommensurate potential [13| . After counting 
on the locations of Anderson plateaus, we find that the 
height of Anderson plateaus is decided by a' with values 
a', 1-a', 2a', l-2a', 2(l-a'), l-2(l-a'), 4a', l-4a', 
if the values are in the range of (0,1). The Anderson 
plateaus with hi — a', 1 — a', 2a', 1 — 2a' are the most 
important as they are wide enough and easy to appear 
for Vj < 2. So heights of Anderson plateaus are totally 
decided by a', while the width and the existence of the 
Anderson plateau are associated with Vi and the charac- 
teristic density (/J). For system with high characteristic 
density, except of the existence of Anderson plateaus, 
there is a Mott plateau in the trap center characterized 
by m = 1. 

Finite-temperature density distribution and quantum 
criticality.- The finite-temperature density distribution 
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can be calculated by 
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ni (T) = -J2e- E ' 
n=l 



/k B T 



n C B\b}bi\* 



HCB/ 



(3) 



where N s = L\/{L — N)\N\ is the number of states, E n is 
the energy of eigenstate $g CB , and Z = J2n=i e~ En / kBT 
is the canonical partition function. By using the exact 
numerical method in Ref. (l9j , we can calculate the finite- 
temperature properties of the hard-core bosons in incom- 
mensurate lattice very efficiently. In Fig|3l we display 
the local average density profiles for systems with differ- 
ent temperatures. The obvious plateau at zero temper- 
ature tends to vanish with the increase in the temper- 
ature. Consequently, the sharp boundary for supcrfluid 
and Anderson localized phases is destroyed by tempera- 
ture fluctuations. If the temperature is low and located in 
the quantum critical regime, the density distributions for 
different temperatures should fulfil some universal scaling 
laws around the quantum phase transition point. 

Before going to the analysis on the quantum critical- 
ity, we would like to introduce briefly the general theory 
for the quantum criticality. Given the equation of state 
ii = n(T, /x) for a system with density n, temperature T 
and chemical potential /x, then near the zero-temperature 
phase transition point /x = /x c , it was shown (jl. Il5| that 
when the dimensionality of the system is below a critical 
dimension, the following universal relation exists: 



n(/i,T) - n r (/x,T) = T" 



-£n(£5&), ( 4 ) 



with n,-(/x, T) being the regular part of the density, d the 
dimensionality of the system, v the correlation length 
exponent, and z the dynamical exponent. Q(x) is a uni- 
versal function which describes the singular part of the 
density near criticality. If n r (/x, T) in the above equa- 
tion is known for one or more phases (such as vacuum 
and Mott insulator), Eq.Q would be very useful for de- 
termining the quantum phase transition point. After 
scaling the density with the T term, we get A(n,T) — 

T~^~ 1+ ~^(n{fi,T) - n r (n,T)). If we plot A(fi,T) ver- 
sus (i for a system at different temperatures, then curves 
with different temperatures will intersect at the same 
point (i = fi c . With this scaling method, we can de- 
tect critical points of phase transitions at zero tempera- 
ture from none-zero-temperature density profiles. Once 
/x c is determined, we can scale the chemical potential 

~ -L ~ 

fx = (fj, — [ic)/Tvz. By plotting A(/j,,T) versus ji for 

a system at different temperatures, all curves collapse 

into a single one (the universal function VL{x)). It was 

shown that d = 1, z = 2 and v = 1/2 for both the 

vacuum-superfluid and superfluid-Mott phase transition 

of hard-core bosons 



15|. 



For hard-core bosons in the ID incommensurate lat- 
tice within a harmonic trap, there may exist vacuum- 
superfluid, superfluid-BG, and even the superfluid-Mott 




FIG. 3: (Color online) Local average density profiles for sys- 
tems with different temperatures. Insert: enlargement of the 
Anderson-plateau area. The system is with 1000 sites, 200 
bosons, Vi = 1.3, a = (\/5-l)/2 and V H = 2.5 x 10" 5 . 
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FIG. 4: (Color online) (a,b): Scaled local average density 
profiles vs fi(i) at different temperatures. (c,d,e): Scaled lo- 
cal average density profiles vs scaled chemical potential at 
different temperatures for jU c i(c), /-t C 2(d), /-i C 3(e). (f): The 
local average density profile of the system at zero tempera- 
ture together with the critical point i c got from the scaling 
method. The system is with 1000 sites, 200 bosons, Vi = 1.3, 
a = (VE- l)/2 and Vh = 2.5 x 10~ 5 . 



insulator phase transition as the position i changing from 
sides to the center. If we replace the chemical poten- 
tial fi in Eq.((4|) by fx(i) = fx — Vh(i), where V/j(i) is 
the harmonic trap potential and /x is determined by 
^2iTi(iJ,(i),T) = N, then /x(j) will drive the transition 
from one phase to the other phase. Now we explore the 
quantum criticality of the disordered Bose system by us- 
ing the local average density distribution. We find that 
there exist similar scaling laws described by Eq.(j4]) but 
with the density distribution in Eq.(j3]) being replaced 
by the local average density distribution. To give con- 
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FIG. 5: The right half of the phase diagram at zero temper- 
ature in Vj-position plane of the system with 1000 sites, 600 
bosons, a = (y/E - l)/2 and Vh = 2.5 x 10" 5 . 



crete examples, we consider a system with 200 hard-core 
bosons, Vi = 1.3, a= (\/5- l)/2 and V H = 2.5 x 10~ 5 , 
for which there is no Mott insulating plateau in the 
local average density profile as the characteristic den- 
sity p = 1 < p c . In Fig we plot scaled densities 
f%i(T) /T 1 / 2 for different temperatures around the edge 
of the density profile, where fii(T) = n(pi,T). Differ- 
ent curves intersect at the point p c i(i) = —2.293, which 
is just the zero-temperature vacuum-superfluid transi- 
tion point and occurs at sites i c \ = 788.2 and 212.8 via 
fx(i) = p — Vh (i) ■ 

From Figj3]we know that Anderson plateaus appear at 
fii = a 1 = 0.38197. By taking the regular part n r = a' , 
in Figfjb we plot scaled densities (nj(T) — a') /T 1 / 2 versus 
p(i) around the regime of Anderson plateau for different 
temperatures. It is clear that the different curves inter- 
sect at two points, which indicates that the superfluid- 
BG phase transitions occur at p C 2{i) = —1.472 and 
Hc3(i) = — 0.258or equivalently at i c2 = 724 and 277, 
and at i c s = 537.6 and 463.4. To display the universal 
scaling functions, we plot scaled local average densities at 
different temperatures against scaled chemical potentials 
around critical points p c \, p C 2, p C 3 in Fig|4j:, d, e. We 
observe that curves for different temperatures indeed col- 
lapse onto a single curve after scaling, except that there 
are oscillations on the curves caused by the incommen- 
surate lattice. In Fig|4f we show the zero-temperature 
local average density profile in contrast to critical points 
determined by the scaling method and marked by dash 
lines. 

Through the analysis of quantum criticality, we can 
map out the zero temperature phase diagram from the 
nonzero-temperature local average density distributions. 
As an example, we show the phase diagram of the sys- 
tem with 1000 sites, 600 bosons, a — (y5 — l)/2 and 
Vh = 2.5 x 10~ 5 in the V/-position plane, which has the 
characteristic density p = 3 > p c . As shown in FigJSJ 
without the incommensurate potential the system has a 



structure with the Mott insulator in the center surround- 
ing by the superfluid phase on two sides. As Vi increases, 
Anderson plateaus appear in the local average density 
profiles and the particles in them are in BG phase. With 
further increase in Vj , more plateaus appear and become 
wider. Furthermore, the size of Mott insulator decreases 
but the system size grows slowly as Vj increases. 

In summary, using the Bose-Fermi mapping, we study 
the universal properties and quantum criticality of one- 
dimensional hard-core bosons on incommensurate opti- 
cal lattices within harmonic traps. By calculating finite- 
temperature density distributions exactly, we unveil uni- 
versal scaling relations of local average density distribu- 
tions in the quantum critical regime. The zero temper- 
ature phase diagram is then determined from nonzero 
temperature local average density profiles with help of 
the analysis of quantum criticality. 
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